Exact wave functions for concentric two-electron systems 
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We show that the exact solution of the Schrodinger equation for two electrons confined to two 
distinct concentric rings or spheres can be found in closed form for particular sets of the ring 
or sphere radii. In the case of two concentric rings, we report exact polynomial and irrational 
solutions. The same methodology is applied to the case of two concentric spheres for which we 
report exact polynomial solutions for the ground state and the excited states of S symmetry. For 
these concentric systems, we show that the exact wave function does not contain terms proportional 
to the interelectronic distance due to the spatial separation of the electrons. 
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I. INTRODUCTION 



In the past decade, much effort has been devoted to the fabrication of GaAs/AlGaAs concentric double quantum 
rings using droplet-epitaxial techniques [lj-ll] in order to study the Aharonov-Bohm effect [1] and the influence of the 
Coulomb interaction on the magnetic properties More recently, accurate theoretical calculations, using local spin 
density functional theory exact diagonalization Q and quantum Monte Carlo [7j techniques, have shown that 
Wigner molecules [1] are formed at low density. 

Multiple-shell fullerenes or "buckyonions" 0, [l(| are another kind of concentric material which have attracted 
research attention for many years. Since the 1985 discovery of fullerene [ll|, experimental search for such carbon- 
based materials have been very fruitful, leading to the discovery of carbon nanotubes, boron buckyballs, etc. The 
morphology of buckyonions and the interaction energy between shells have been investigated using first principles 
calculations [IH, [l3[ and continuum models 0, [l5[ . 

A minimalist model for these two concentric systems consists of electrons confined to two concentric rings or spheres. 
In Ref. [la ] , we reported a comprehensive numerical study of the singlet ground state of two electrons on concentric 
spheres with different radii. We analyzed the strengths and weaknesses of several electronic structure models, ranging 
from the mean-field approximation (restricted and unrestricted Hartree-Fock solutions) to configuration interaction 
expansion, leading to near-exact wave functions and energies. Berry and collaborators also considered this model to 
simulate doubly-excited states of helium [ljj and the rovibrational spectra of the water molecule in both the ground 
[l8j and excited states [l9[ . More recently, the model has been applied to quantum- mechanical calculations of large- 
amplitude light atom dynamics in polyatomic hydrides [20l l2l| , confirming again the broad applicability of theoretical 
models of concentric electrons. 

In a recent series of papers [22U24]] . we have found exact solutions for several new quasi-exactly solvable models 
[25l | . These are valuable for the testing and dev elop ment of approximate methods, such as exchange and correlation 
functionals within density-functional theory 26 -28]. In particular, we have shown that one can solve the Schrodinger 
equation for two electrons confined to a ring [24| or to the surface of a sphere [22|, [H| for particular values of the 
ring or sphere radius. Subsequently, Guo et al. reported some exact solutions for two vertically coupled quantum 
rings using the same technique 29] but their analysis is incomplete. Although the present study focuses on different 
systems, our analysis can be applied to find additional exact solutions for their system. 

In this Article, we extend our analysis to the case where two electrons are located on concentric rings or spheres 
with different radii. We focus our attention on the rotationally-invariant (i.e. S) states. Although our principal 
concerns are the T> = 1 (electrons on rings) and T> = 2 (electrons on spheres) cases, we consider the general case in 
which two electrons are confined to two distinct 2?-spheres of radius R\ and R2 30]. For R\ 7^ i?2, the singlet and 
triplet states are degenerate and one can easily switch from singlets to triplets (or vice versa) by antisymmetrizing 
(or symmetrizing) the spatial wave function [la] . For T> = 1, this degeneracy is observed even when R\ = R2 because 
the strength of the Coulomb singularity at u = creates a node there [13] ■ 



II. EXACT SOLUTIONS 



The electronic Hamiltonian on the 2?-sphere is 



a~3±3 + i, (1) 

2 u 



where u = ^fR\ + R\ - 2R X R 2 cosw is the distance between the two electrons, V 2 is the Laplacian in hyperspherical 
coordinates [3lf , and u £ [0, n] is the interelectronic angle and satisfies 

fcos(0i-0 2 ), £> = 1, ,„s 

cosw = < (2) 
I cos 0\ cos 02 + sin 0i sin0 2 cos(<pi — 4>2), T> = 2, 

where <pi € [0, 27r) and 0i £ [0, n] are the usual azimuthal and polar angles of the ith electron. For notational simplicity, 
we introduce 

\ = h = h x +^ a ± =nl = ( Rl ±R 2 f, (3) 

where 71 is an "effective" radius, and 1Z- and 1Z+ are respectively the lowest and highest possible values of u. 
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For states of zero angular momentum (S states), the wave function depends only on u and the Hamiltonian becomes 



H = — - 

4cr 



(a- + (T 4 



■^ + r (2P - 1)tt -MMa)-^! 

du z u u A 



d_ 

du 



1 

u 



The resulting Schrodinger equation is Fuchsian [331 ] with singularities at —1Z+, 0, 1Z- and 1Z+: 



*"(u) 



2?/2 



D/2 



1 



D/2 



D/2 



it + 1Z+ u + 1Z- u u — 1Z- u — lZ-\ 



V(u 



(u 



(4) 



(5) 



where V(u) = 47Z 2 u(l — Eu) is a Van Vleck polynomial [Hj]. The Fuchsian equation is a generalization of the Heun 
equation [34] (which appears when R\ = i£a) for higher number of singularities f22l424|. 
Knowing that Eq. ((SJ is Fuchsian, we seek solutions of the form 



2 \ 1/2 



2 \ 



72 2 



where (a, 6) can take the values (0, 0), (2 - 2?, 0), (0, 2-2?) and (2 - 2?, 2 - 2?). The functions 

n 



(6) 



(7) 



fe=0 



are called nth-degree Stieltjes polynomials [3J]. The index m is the number of nodes between 1Z- and 72.+ . The 
different possible values for the set (a, 6) are determined by ensuring that the substitution of ^ into ([3]) conserves 
the Fuchsian nature of the differential equation and the Van Vleck nature of the polynomial V(u). After making this 

i{a,b) 



with 



1 

47 



p{u) fS^u) + {aib) dS%£'(u) + S%£>(u) = E (a,b) S (a£) {u)t 

du 2 du u 



-r(a,6)/ 



,(o,&) , 



P(u) = u 2 + — ^t-( a _ +(T+ ), 

(u) = [(27? - 1) + 2(a + b)} u - <P~ + °+) + ^ 

u u 6 

(0, (0,6) = (0,0), 

^Kb) = S + J 23(2? - 2)/(4d), (a, 6) = (2 - 2?, 0) or (a, 6) = (0,2- 2?), 
((2?-2)/a, (a, 6) = (2 -2?, 2 -2?). 

It can be shown that the coefficients satisfy the five-term recurrence relation 



(8) 

(9) 
(10) 

(11) 



Sfc+4 = 



(fc + X>)(<7_ + 0+) +2(acr_ +6(7+) 



(7_<7+ k + 4 



with the starting values 

so = 1, 



4a 



4cr 
fc + 2 



*i = 0, 



-Sfc+2 



Sfc+1 + 



fc[fc + 22? + 2( a + b-l)] _ M 
4cr 



4a 



S3 = 



3(7- (7+ 



ST.- 



(12) 



(13) 



One can note the similarity of the three-term recurrence relation obtained for R\ — R2 [22j-|24| and the term in square 
brackets in (|12l) . The condition Sj = shows that the exact wave functions do not contain any term proportional to 
the interelectronic distance u. Indeed, due to the spatial separation of the electrons, there is no need to satisfy the 
electron-electron Kato cusp condition (35j . 
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TABLE I. Solutions for the ground and excited states of two electrons on concentric rings (T> = 1). 



m 


n 


(a, 6) 


Ri 


i?2 


E 





4 


(0, 0) 


6.56723 


15.0297 


0.0552268 





5 


(0,0) 


9.69674 


14.7302 


0.0476376 





6 


(0,0) 


14.5140 


18.9048 


0.0339531 





7 


(0, 0) 


20.4922 


24.6415 


0.0246731 





8 


(0,0) 


27.5378 


31.5687 


0.0185769 


1 


5 


(1,0) 


17.2072 


47.7424 


0.0171724 


1 


6 


(1,0) 


20.3910 


37.9387 


0.0189862 


1 


7 


(1,0) 


26.3702 


40.6719 


0.0163405 


1 


8 


(1,0) 


33.8677 


46.8166 


0.0134467 


1 


4 


(0,1) 


8.10699 


19.293 


0.0559434 


1 


5 


(0,1) 


10.9336 


16.1063 


0.0549899 


1 


6 


(0,1) 


16.0422 


20.4823 


0.0384000 


1 


7 


(0,1) 


22.3238 


26.4996 


0.0274453 


1 


8 


(0,1) 


29.6705 


33.7181 


0.0204070 


2 


7 


(0, o) 


12.1683 


17.5157 


0.0613308 


2 


8 


(0, o) 


17.5501 


22.0442 


0.0424363 


2 


5 


(1,1) 


19.9728 


57.0944 


0.0172333 



l, 


u 


(~\ 1 ~1 
U, L ) 


99 9 c i c in 




n D91 ftd^R 


2 


7 


(1,1) 


28.4414 


42.9590 


0.0180032 


2 


7 


(1,1) 


36.2331 


49.2980 


0.0146648 


3 


7 


(1,0) 


24.9101 


81.0762 


0.0141096 


3 


8 


(1,0) 


24.1637 


43.1663 


0.0227745 


4 


8 


(1,1) 


26.1359 


46.2032 


0.0241549 



To define the nth-degree polynomial Sn^n completely, we need to find the values of the inner and outer radii R\ 
and i?2, the exact energy E, and the coefficient S2- These are obtained by solving s„+i = s„+2 = s n +3 = Sn+4 = 0, 
which arc functions of these four unknown quantities. For s n+1 = s n+2 = s n+3 = 0, the condition 

E (a,b) = n[n + 2V + 2(a + b-l)] ^ 

ensures that s„+4 = and allows to determine the exact energy knowing the degree of the polynomials and the 
effective radius. 



III. CONCENTRIC RINGS 



For V = 1, i.e. two electrons confined to two distinct rings of radius R\ and i?2, we have four different sets of 
solutions corresponding to the values (a, b) — (0,0), (1,0), (0, 1) and (1, 1). They correspond to the wave functions 

= *£#(u) = cos !s&°>(u), (15) 

= sin |s£#(u), = sinu S^(u). (16) 

As shown in Table |U the family (0, 0) produces ground states and excited states wave functions, while the other 
families only produces excited states with at least one node (m > 1) for (1,0) and (0, 1) and two nodes (m > 2) for 
(1,1). 

The first nodeless solution is 



S^(u) 1 - ————V? - -^-u* - 16807 u * (17) 
\o W-l 1521 « n8638 " 240597864"' (17j 
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TABLE II. Solutions for the ground and excited states of two electrons on concentric spheres (T> = 2). 



m n Ri R.2 E 

5 13.4767 27.3109 0.0299540 

6 18.0008 27.7797 0.0262917 

7 24.2569 32.8421 0.0206850 

8 31.7451 39.7970 0.0162370 

1 7 17.0759 34.7323 0.0335356 
1 8 21.3402 31.6212 0.0319596 



for 



i?i = y ^3(13-^78), i? 2 -y\/3(13 + V78), (18) 

and the energy of E = 28/507 = 0.0552268. 

The first solution which has a sole node at w — is produced by the family (0,1) for n — 4 and reads 

5(0.1) () i250_ _ 15625 _ 390625 

^ ( > ~ 58653 9853704 M 12231731232 ' 1 9j 

for 



#i = ^399(7- 2^6), i? 2 = — ^399(7 + 2%/6), (20) 



and the energy of E = 625/11172 = 0.0559434. 



IV. CONCENTRIC SPHERES 



In the case of two electrons on concentric spheres (V = 2), there is only one family of solutions and it is associated 
with (a, b) = (0, 0). The first few exact solutions are reported in Table ITTT The first solution is obtained for n = 5 and 
is a ground state wave function 

S 5,o\ u ) = 1 - 0.0161923 u 2 - 0.0012233 u 3 - 0.000033429 u 4 ~ 3. 14705 x 10~ 7 u 5 , (21) 

associated with the energy E = 0.029954 and the radii i?i = 13.4767 and R 2 = 27.3109. 

Wave function terms proportional to u are difficult to model using configuration interaction (CI) methods, and it is 
common to use explicitly correlated methods [36j to improve the convergence of such calculations [37] . This was well 
illustrated in our i?i = R 2 calculations [Hj], where we observed slow convergence of the CI expansion with respect to 
basis set size. In contrast, the absence of a linear term in the wave function pip explains the rapid convergence of 
our Ri 7^ i? 2 calculations fig . 



V. CONCLUSION 



In this Article, we have shown that the Schrodinger equation for two electrons confined to two distinct concentric 
P-spheres can be solved exactly for particular sets of the P-sphere radii. In the case of two concentric rings (D = 1), 
we have reported the four sets of exact wave functions. We have observed that the exact wave functions do not 
contain any term proportional to the interelectronic distance u. In this case of two concentric spheres (T> = 2), we 
have reported exact ground- and excited-state wave functions. 

We have not reported any exact solutions for T> > 3 in the present manuscript, but these are not difficult to derive. 
However, we note that only the (a, b) = (0,0) family is physically significant because, for (a, b) =t (0,0), the exact 
solutions of © are divergent at u = or u = 1Z + . 

A natural extension of the system studied in this Brief Report would be to consider two rings of different radii 
vertically separated by a distance d. We believe this system is also quasi-exactly solvable, but the derivation would 
be more complicated because of the larger number of singularities in the Fuchsian equation. For this reason, we chose 
not to consider this system in the present study. 
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